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LIMITS OF VALIDITY OF RADIATIONLESS DECAY OF HIGH-DENSITY PLASMAS 


by Norbert Stankiewicz 
Lewis Research Center 

SUMMARY 

The re combinative decay of collision-dominated plasmas is investigated theoretically, 

and the limit of validity of such radiationless models is determined. The density range 

considered conforms to those plasma densities usually encountered in the analysis of a 

magnetohydrodynamic generator. The results indicate that neglecting radiative processes 

may not always be justified at these densities. It is concluded that the radiationless model 

15 

is valid for argon when the electron density is greater than 6. 8x10 per cubic centimeter 

15 

and for cesium when it is greater than 1 . 1x10 per cubic centimeter. 

INTRODUCTION 

One of the problems in plasma physics, and in particular, magnetohydrodynamic 
(MHD) power generation is the determination of the range of validity of the Saha equation. 
This equation gives the fractional ionization of a gas in equilibrium as a function of neu- 
tral density, temperature, and ionization potential. It is common practice in MHD gener- 
ator analysis (e. g. , ref. 1) to apply this relation to the nonequilibrium situation in which 
the plasma electrons are at a higher temperature than the heavier species. This temper- 
ature difference is caused by the electron heating phenomenon produced in a plasma flow- 
ing through a magnetic field (ref. 2). The electron density is normally calculated from 
the Saha equation at this higher electron temperature. 

Previous work on this problem includes the hydrogen calculations of UPyanov (ref. 3) 
using Born approximation cross sections; the model atom (based on cesium) calculations 
of BenDaniel and Tamor (ref. 4) using cross sections based on the Townsend a coeffi- 
cient; and cesium and potassium (in helium and argon atmospheres) data of Hiramoto 
(ref. 5) based on the J. J. Thompson theory of recombination (ref. 6); and the work of 
Griem (ref. 7) using semiclassically derived cross sections. 

This study investigates the Saha assumption for a decaying (i. e. , recombining) 
plasma, whose electrons are initially at a higher temperature than the ions and neutrals. 


The initial conditions are assumed to be such that the dominant mode of plasma relaxa- 
tion is by three-body recombination in which the third body is another electron. In this 
regime, the bound states are in equilibrium with the free electrons and the Saha equation 
therefore holds. The conclusion that optically dense plasmas (collisionally dominated) 
obey the Saha relation at the electron temperature has been previously reached in the the- 
oretical discussions of reference 8. Experimentally, this behavior has been observed in 
reference 9 for a highly ionized hydrogen plasma Likewise, the argon data presented in 
reference 10 indicate that for the higher electron densities the decay follows the Saha re- 
lation. 

As the decay continues and the electron density decreases, collisions with the neutral 
atoms are not frequent enough to maintain this equilibrium state. Radiative processes 
then become competitive with collisional effects in determining the bound level populations 
An estimate of this threshold point is made based on the assumption that the ground state 
of the atom is the first level to fall out of equilibrium with the free electrons. The 
Gryzinski expression for the inelastic cross sections (ref. 11), is used and the threshold 
point is calculated for argon and cesium. 


ANALYSIS 
Basic Premises 

As in references 3, 4, and 7, electron-electron interactions are assumed to occur 

so rapidly that a Maxwellian velocity distribution of electrons is assured. Accordingly, 

3 3 

the number of electrons in a six-dimensional phase cell dq e dv g is given by 


3 3 / „ \3/2 9 

£ e dq e dv e = "e (i^f) ex P(-"V 2kT e> <*> 


3 

where f is the electron distribution function; T , the electron temperature; dq , the 
6 ^ ^ 

infinitesimal volume in coordinate space; and dv g , the infinitesimal volume in velocity 
space. (All symbols are defined in appendix A. ) The distribution function in equation(l) 
is normalized to the electron number density n g by integrating over all velocity space. 

Hence, 

"e -/< [ e 4 < 2 > 


2 



The local electron thermal energy u g is given by a similar integration: 

"e-/| mv l f e4-| n . Kr e <3 > 

It, 

The number density of neutrals in the j n excited state is denoted as n., and each is 
regarded as an independent species of the plasma. Since the total number density of nu- 
clei N is a constant, the number of independent equations describing the plasma is 
reduced by one because of the normalizing equation: 


E n j + "i - N =Z n ) + “e < 4 > 

j j 

The last step is a consequence of the assumed charge neutrality of the system. (Only 
singly charged ions are considered. ) 

If the neutral levels happen to be in equilibrium with the free electrons, each excited 
state will obey a Saha relation of the form 


s 2 
n j ~ n e 


,3/2 


2<Uj \27rmkT e/ 


co. exp(E./kT e ) 


( 5 ) 


The superscript s is used to denote Saha equilibrium; to. is the ion degeneracy and 

1 th 

co. and E . are the degeneracy and the ionization energy of the j state. 


General Formulation of Problem 


The following processes are considered to be most likely for the plasma densities of 
interest in this study: 

i-U 

(1) Three-body (electron-electron-ion) recombination to the ] 1 state and its inverse 
two-body (electron-neutral) ionization from the j n state: 


e+e + i^e+n. 

(2) Two-body (electron- ion) recombination with photoemission and its inverse, photo- 
ionization: 


e + i 


:n j + 


hv 


3 


(3) Superelastic and inelastic collisions: 


e + n^ ~ e + n k j =£ k 

(4) Radiative deexcitation and photoexcitation: 

n k r n j + hv k > i 

(5) Recoil energy loss by the electrons in elastic collisions with the ions and neutrals: 

ei+ (n)~ e 2 E i > E 2 

For these processes, the coupled rate equations describing the decaying plasma can be 
written as 



The superscripts c and r refer to collisional and radiative processes, respectively, 
while the subscripts IR, j~k, and N refer to ionization recombination processes, tran- 
sitions between levels, and elastic recoil energy losses that the electrons suffer during 
collisions with the ions and neutrals. The differential equation describing the tempera- 
ture decay is obtained by differentiating equation (3) and substituting equations (6a) and 
(6c): 


4 


( 7 ) 



j,k 


The following forms of the collisional terms in the rate equations (6) are derived in 
appendix B: 


(^) = n e E ( n j ' n f) <vo I< E j “ E *» 

($) = "ef”! ' n l) <TO I (E J " E * )> 

\ / IR N 7 

^ n j ex Pt (E k “ E j)/ kT e| (va ex( k ~ «> 


k<j 


+ nA n k exp[(E j - E k )/kT e ] - n. 


k>j 


©L^X ( n ? - n Kl mv2vCT i (E 3 - e *>> 


= n °)' j Y J nk 6XP[(Ej " Ek)/kTel ' n )t mv2va ex ( ^ k) ) 

i,k i k>i 


i k>j 


From reference 1, the recoil energy losses can be written as 


( 8 ) 


(9) 


KvCT e x(j -k)> (10) 


( 11 ) 


( 12 ) 


m 


— 1 = - — n e k (T e - T 0 )(^ ei + * en ) 


, 9t /N M 


(13) 


5 



where Tq is the temperature of the ions and neutrals (taken to be equal), v ^ and y gn 
are the electron- ion and electron- neutral collision frequency for momentum transfer, m 
is the electron mass, and M is the mass of the neutral (or ion). 

Finally, from equation (4), it is clear that 


dn 


e 


dt 




(14) 


Radiationless Case 

For the radiationless case, consideration of equations (6), (8), (10), and (15) shows 

that 



(15) 


This equation implies that the inelastic-superelastic collision rates vanish by detailed 
balancing. Equation (10) then yields 


— = — exp[(E - Ej/krj E > E (16) 

n. co. ] k e J k 

J J 

which indicates that under these conditions, equilibrium exists among the bound states at 
the free-electron temperature. Consequently, the Saha relation in equation (5) holds for 
each bound state. 

If equation (5) is summed and inserted into the normalizing condition, equation (4), 
the result gives the usual form of the Saha equation: 


n 


2 h L 1 

20^(2 7rmkT 


1 3/2 


exp(E Q /kT e )^qj j exp[(E Q - E-J/kTJ = N - n g (17) 


For moderate temperatures, the summed term (the partition function) is equal to the 
ground- state degeneracy oj q and therefore equation (17) may be written as 


6 



(18) 


n 


N - n 


J^(27raikT e ) 3/2 exp(-E c /M' e ) 

co h 3 
o 


For the radiationless case then, the electron density is only a function of the electron 
temperature, and a complete description of the decaying plasma is formulated by solving 
equation (7). Because of the equilibrium conditions given by equations (16) and (5), equa- 
tion (7) reduces to 




and substituting equation (13) gives 



~( T*-T n )(i/ ei +0 


M 


en' 


(19) 


( 20 ) 


Although the results given in equations (17) and (20) are well known, the proof pre- 
sented here is useful for determining the extent that a competing process, such as the 
loss of energy through radiation, causes a departure from detailed balancing, as will be 
shown in the next section. 


Conditions of Validity of Radiationless Case 

As the temperature decreases, the number of electrons capable of maintaining the 
superelastic and inelastic processes in detailed balancing also decreases. The ground 
state will be the first level to depart from Saha equilibrium primarily because of its 
larger energy separation from the continuum. Therefore, in order to estimate the equi- 
librium threshold, all states above the ground state are assumed to remain in equilib- 
rium; that is, 


n j = j * 0 (21) 

All emitted radiation arising from transitions to the ground state is also assumed lost, 
that is, the plasma is optically thin. 


7 



Since the total gain in neutral population in all levels is equal to the rate that free 
electrons leave the continuum, the following is true: 



j 


Then, from considering equations (7) and (15), it is concluded that 



(23) 


It is no longer valid to argue that equation (23) holds because of detailed balancing, since 
by hypothesis some inverse processes are missing. In particular, for a transparent 

/dn.\ r 

plasma the term [ — 1 ] contains no reabsorption term from the ground state. For the 


states above ground, since they are in Saha equilibrium and therefore collisionally domi- 
nated, dismissing radiative processes between them as being improbable is consistent. 
The collisional terms vanish between states that are in equilibrium, as seen from equa- 
tion (10), and equation (23) then becomes 



(24) 


The first term in equation (24) is given by equation (10) and is equal to 



= n 


o~k 



- n k exp[(E Q - E k )/kT] - n Q f <va ex (o-k)> 
k 


(25) 


Since the states above the ground level are assumed to be in equilibrium, use of 
equation (6) yields 


8 



( 26 ) 


CO 


3 / \3/2 

~ < “**• - V*. ■ "e sfey “o - < 


Hence, equation (25) can be written as 


L 


^7 o=k =n °( n °' n °) Z <v<w ° * k)> 


(27) 


where V* (va (° — k)> is the collision coefficient for inelastic collisions from the 

GX 

k 

ground state to all excited states. 

The radiative deexcitation term in equation (24) is proportional to the population of 

i. L 

the k in state; the constant of proportionality is the spontaneous emission probability to 
the ground state A(k — o). Hence, 



9 Z 

0t 


o— k 


= ^ njA(k - 


o) 


(28) 


k>o 


k>o 


Using equation (26) (with o> 0 = 1) will change equation (28) to 

E f^ 2 ) = "‘Z“k ex rt‘< E o - E k )/kT e ]A(k - o) 

v t / o— k k>o 


(29) 


Substituting equations (29) and (27) into equation (24) yields the ratio of the nonequil- 
ibrium ground-state density to the equivalent Saha density of the ground state: 


Zu i + .^° 

s 

n 


2 ^k ex PH E 0 ■ E k^ kT el A ^ k - °) 


n e E < v e CT ex(° ~ k » 
k>o 


(30) 


When equation (B23) is substituted into equation (30), equation (30) becomes 


9 



(31) 


ofe exp[-(E Q - E k )/kT e ]A(k - o) 

n o n e E exp H E 0 “ E k)/ kT eK v eW k * °» 

k>o t 

The second term of equation (31) is the ratio of the total radiative deexcitation rate 
to the total collisional deexcitation rate into the ground state. If the collisional rates 
greatly predominate, equation (31) shows that the ground state is in equilibrium, as ex- 
pected. 

For the gases and densities of interest, the contribution to the series in equation (3 1) 
from terms higher than the first is smaller by at least an order of magnitude. This 
change in contribution is due to the reduced population of higher levels, the smaller cross 
sections, and the smaller radiative transition coefficients. 

Equation (31) then reduces to 


n 

= 1 


z 

k>o 


n„ 

-2=l + . 
nf 


A(1 - o) 

n e< va dex^ ~ °» 


(32) 


However, in order to use the summed rate coefficients as calculated in reference 9, 
equation (32) is finally written as 


= j | A(1 — o) 

n o n e £ <TO dex< k “ °» 
k>o 


The error in passing from equation (32) to equation (33) is about 10 percent. The 
effect of increasing optical thickness, however, is in the direction of reducing this error. 

The conductivity of a plasma is of primary interest in MHD generators. Since the 
conductivity a is directly proportional to the electron number density, the fractional de- 
crease in conductivity from its Saha value can be written as 


a 


s 


a 


- g 

s 



( 34 ) 


In the region where deviations from the Saha values become important, it is expected 
that most of the neutrals are in their ground state. Equation (34) then becomes 


10 



( 35 ) 



Substituting equation (33) and defining F? as the fraction of ionization calculated 
under Saha conditions allows equation (35) to be written as 


A(1 - o) t 1 ~ F i> 

n e 2 <VOr dex< k ~ °» Fi 
k>o 


(36) 


Equations (33) and (36) form the basis of the next section in establishing a criterion 
for equilibrium. 


RESULTS AND DISCUSSION 

Equation (33) may now be applied to determine the threshold of validity of the colli- 
sionally dominated plasma for those conditions (pressure and temperature) that are re- 
quired in MHD generators. Reference 1 indicates that argon seeded with cesium is a 
useful working fluid for such a generator. The remainder of this report will therefore be 
concerned with determining the equilibrium threshold for these two gases. 

From reference 12 the spontaneous transition coefficient A(1 — 0) for cesium is 
7 

3x10 per second. Dugan (ref. 13) calculated the superelastic collision frequency for 

15 

cesium for an electron density of 10 per cubic centimeter using Gryzinski cross sec- 
tions (ref. 11). Tables II and HI of reference 13 give these data, from which. the term 

^ ^ va dex^ k °)) ma Y be obtained by dividing by 10 15 . The resulting values are 
k>o 


r 


9. 54x10” 8 cu cm/sec; T = 1000° K 


2 < VCT dex< k ^ °» = 
k>0 


< 

7. 64X10” 8 cu cm/sec; T = 5000° K 


Since A(1 — o) is not a function of temperature, the ratio 


(37) 


11 



A(1 - o) 

£ < TO de X < k -°» 
k>o 

does not change drastically with temperature. The threshold being sought is expected to 
be nearer the 5000° K data point. Equation (33) then gives 


n o- n o 3. 93X10 14 


n 


n 


(38) 


The criterion for equilibrium will be taken arbitrarily as a deviation (n Q - n®/n® 

from the ground state density of 1/e or less. Equation (38) therefore shows that for elec- 

15 

tron densities above 1 . 1x10 per cubic centimeter, the cesium ground state is adequately 
given by the Saha expression. 

For argon, Dugan (personal communication) estimated the resonance transition co- 

O 

efficient A(1 — o) by extrapolation from helium and neon data to be 5x10 per second. 

The computer program for superelastic collisions used by Dugan in the cesium calcula- 
tions of reference 13 was adapted for argon and the following information was obtained: 


£ < w dex <k -°» 
k>o 


2X10" 7 cu cm/sec; T = 10 000° K 
< 5. 6x10” 7 cu cm/sec; T = 5000° K 
1. 6x10" ® cu cm/sec; T = 2000° K 


Using the 10 4 0 


K data changes equation 


(33) for argon to 


(39) 


TABLE I. - EQUILIBRIUM LIMITS 


Gas 

Electron number density, electrons/ cm' 


Present 

report 

Griem 
(ref. 7) 

BenDaniel 
and Tamor 
(ref. 4) 

Chen 
experi- 
mental, 
ref. (10) 

Argon 

6. 8xl0 15 

8. OxlO 16 


-3X10 15 


Cesium 

1. IxlO 15 

1. IxlO 16 

~1. 2X10 17 




n o ' n o 2. 5X10 15 


n 


n 


(40) 


Applying the arbitrary deviation of 1/e 
shows that the ground state is essentially in 
equilibrium for a free-electron density 
greater than 6. 8xl0 4 ^. 

Comparison of these results with the 
work of other authors is summarized in 
table I. The entries refer to electron num- 


12 



ber densities below which equilibrium should not be expected. These values are compar- 
able to the electron number densities typically used in MHD generators. However, in 
order to apply these results to actual generator conditions, a criterion based on the 
plasma conductivity rather than on electron number density will be established. Accord- 
ing to reference 1, at an ambient temperature of 700° K and a pressure of 10” ^ atmo- 

o 

sphere, the fraction of seed required was less than 10 . Typical neutral number den- 

18 16 

sities are 10 argon and 10 cesium atoms per cubic centimeter. The electron tem- 
peratures corresponding to the electron densities shown in table I are therefore 9250° K 
for argon and 3000° K for cesium (ref. 14). Since electron temperatures higher than 
3200° K are unnecessary (ref. 15) in an actual generator, it can be concluded that the 
argon carrier gas will not contribute to the conductivity. 

If the criterion for equilibrium is arbitarily taken to be a 10 percent deviation of 
conductivity, equation (36) for cesium can be written 


s _ 
a - a 


3. 93X10 14 ^ ” F i) 


< 0 . 1 


n 


F i 


(41) 


When numerator and denominator are multiplied by N, the cesium seed density, 
equation (41) becomes 


3. 93X10 14 ^ " F P 


N 




< 0 . 1 


(42) 


or 


F. -- 3 - 93x11)15 (,/ 1 + 4N - 
2N \V 3.93+10 15 

I £> 

Substituting the value of 10 per cubic centimeter for N gives finally 

F i > 0. 46 



(43) 


(44) 


or 


n > 4. 6X10 15 
e 


( 45 ) 
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ii min in 


'.i: 

Hence, this result differs from the previous one by about a factor of four. Equation (44) 
refers to the fractional ionization of the cesium seed only. For the same conditions, the 
overall fractional ionization of the cesium-argon mixture is 

F. . . = — > 4. 6X10 15 = 4 55x10 -3 (46) 

h N Ar + N Cs 1.01X10 18 

This is a fraction of ionization that is easily met in most MHD designs. However, equa- 
tion (46) is the criterion for a particular choice of gas temperature and operating pres- 
sure. For other conditions, a different fractional ionization requirement will hold. 

The analysis presented in this report still leaves unanswered the question of how the 
required conductivity will be achieved in a generator that relies on the electron heating 
phenomenon, and in which the initial electron densities are below the minimum values 
quoted here. 


CONCLUSIONS 

This study of the limits of radiationless decay of the high-density plasmas argon and 
cesium indicates that the neglect of radiative processes may not always be justified for 
those plasmas usually encountered in magnetohydrodynamic (MHD) generator analysis. 
Comparison of the population of the ground state with the equivalent equilibrium popu- 
lation (Saha) shows that serious departure from equilibrium occurs (1) for argon when 

the free-electron density falls below 6. 8x10 per cubic centimeter, and (2) for cesium 

15 

when the free-electron density falls below 1. 1x10 per cubic centimeter. The result of 
ignoring radiative effects below these limits would be to overestimate the conductivity 
of an MHD generator to an extent that depends on its fractional ionization. This over- 
estimation is especially true of those generators which are designed to employ the elec- 
tron heating phenomenon to achieve the required conductivity, and which necessarily 
operate initially with low electron densities. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 4, 1966 
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APPENDIX A 


A(k — j) 


b 

3 

dq 


3 

dv 


E 

E r E 


k 


E 


o 


e 



h 

k 

M 


m 

N 


n 


SYMBOLS 


spontaneous radiative tran- n^, n^ 

sition coefficient from 
k th to j th state, sec - ^ n 

o 

impact parameter 

p 

infinitesimal volume in co- 
ordinate space T 

T 

infinitesimal volume in o 

velocity space 

energy of free electron 

th t 

ionization energy from j tn 

and k* k states u e 

ionization energy from 

ground state v 

electronic charge; Napier- 
ian base 2.718 X 

fractional ionization v 


total fractional ionization 

electron distribution func- 
tion 

Planck’s constant 

Boltzmann’s constant 

mass of neutral atom or 
ion 

electron mass 

equilibrium number den- 
sity of atoms 

number density 


CT dex< k -j 


CT ex<j- k ) 


<Tj(E. — E* 


neutral number density in 
and k^ k states 

neutral number density in 
ground state 

reaction probability 

temperature 

final equilibrium tempera- 
ture; ion and neutral 
temperature 

time 

average electron energy 
defined in eq. (3) 

velocity 

unit mechanism (see p. 20) 

l/kT e (appendix B) 

collision frequency for 
momentum transfer 

conductivity 

) collision induced deexcita- 

tion cross section from 
k^ k to j^ k state (k > j) 

collision induced excita- 
tion cross section from 
j th to k th state (k > j) 

) ionization cross section 

from j*'* 1 state with emis- 
sion of electron having 
energy E* 
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<7j^(E* — Ej) recombination cross sec- 



tion for capture of elec- 
tron having energy E* 
into j state 

to 

degeneracy 

<> 

average value 

a/ at 

change due to interaction 


processes 

Subscripts: 


dex 

deexcitation 

e 

electron 

ex 

excitation 

I 

ionization 

IR 

ionization- re combination 

i 

ion 

3 

state 

k 

k th state 


jrk 

■fin 

process between 3 n and 
k th state 

N 

elastic collision 

n 

neutral 

o 

ground state 

R 

recombination 

0 

equilibrium 

Superscripts: 


c 

collisional 

r 

radiative 

s 

Saha equilibrium 

*t 

characterizing free elec 
trons 


APPENDIX B 


DERIVATION OF COLLISION INTEGRALS 

Throughout appendix B the parameter X = l/kT g will be used for ease of writing. 
The superscripts *, f will be used to characterize the electrons involved in the colli- 
sions. 


Two-Body Ionization - Three-Body Recombination 

The energy diagram (fig. 1) depicts the energy levels of a typical ionization- 
recombination process. 

The change in electron density in a time dt due to ioniza- 

jth 


Ionization 


tion from the j state is 


Recombination 


_| Ej | 


Ionization 


-E* 

-0 


— ) dt= /i 

v 9t /jl J 


f e dv e nj dq n P(E j - E*)dE* 


(Bl) 


Recombination 

Ej 


Figure 1. - Schematic diagram de- 
picting energy levels of a typical 
ionization-recombination process. 


in which the argument of f g is the energy E of the ionizing 
electron; E* is the final energy of the emitted electron; and 
P(Ej — E*)dE* is the probability of emitting into the E* free 
state (see eq. (1919), ref. 16). The neutral particle must find 


3 3 

itself in a volume element dq n equal to dq g if such a collision is to be effected; hence, 


3 3 lfiv^ 

dq„ dv = b db E dE dt 


(B2) 


m 


The ionization cross section for emission into an energy E* is then 


CTj(E • - E*) = 2tt J P(E. - E*)b db 


(B3) 


With these substitutions and with the explicit form of f g from equation (1), equa- 
tion (Bl) becomes 
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'dxi 

I 

, at 


\ 3 /2 

l) = 8nmn(-±-) f n.e~ AE ct t (E. — E*)dE* E dE 

L e \2nm ) Si 1 J 


(B4) 


where the integral sign represents a double integration over E* from zero to E - 
and over E from E. to infinity. 

] fV, 

The change in electron density due to recombination to the j state (in which another 
electron is the third body that carries away the excess energy) is given by 


( 3n 


at 


dt 


'jR 


/ t 3f * 3* 3* 3 * 

f I dv e f e dq e dv e n i dq i P < E ~ E j> 


(B5) 


JL, J- f 

where the argument of f"T is the energy E of the electron to be captured into the j 

+ e 4- 

state; f^ has the energy argument E ' and represents the electron that will carry away 

the excess energy of the reaction; and P(E* — E.) is the probability of capture into the 

th J 

j cn state. 

Microscopically, the ionization and recombination processes of equations (Bl) and 
(B5) are related by the energy equation 


E = E^ + E* + E- 


(B6) 


O 

From equation (B5) it is noted that the volume element of the ion dq. must be equal 

3 t 1 

to dq^ if a collision is to be realized. Hence, transforming variables as in equation 

(B2) gives equation (B5) as 


(^) - - Afe 
WiR J 


P(E* - Ej)(16vr 2 m) 2 b* db* b* db^ dt E* dE* E^ dE^ 


(B7) 


Defining, as in equation (1921) of reference 16, the cross section for recombination 


to be 


°r( e * 


E j) = -y p ^ E * - E j^ b * db * b ^ db ^ dt 


(B8) 


yields equation (B7) as 
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3 /* 

647r 2 m 2 nV — \ f exp[-X(E^ + E*)]a R (E :t 
e \27Tm/ * 


E.)E"f dE 1 " E* dE* (B9) 


Equation (B 6 ) can be used to change variables from E^ to E; and when equations 
(B9) and (B4) are added and summed over j, the net change in electron density due to 
ionization- recombination processes is found: 


K 

Vat 


= 87mm 


^ /2 ynJ 3XP (-XE)a I (E j - E*)dE 


IR 


* E dE 


3/2 XE /* 

- 87 rmn 2 ^-A_^ e 3 yexp(-XE)or R (E* - E.)E* dE*(E - E* - Ej)dE 


(BIO) 


Here the integral sign represents, as mentioned before, a double integration over E and 
E* with their appropriate limits. 

The relation between cr R (E* - Ej) and Oj(Ej - E*) is given in equation (1924) of ref- 
erence 16, and in the notation of this report is, 

Ect t (E. - E*) = E*(E - E* - E.)ct r (E* - E.) (Bll) 

1 h 3 "3 1 1 


nj^TOjfEj - E*)> (B12) 

where the definition of the Saha density in equation (5) was used and 

3/2 

<VCTj(E. - E*)> = 877 m^-^ f exp(-AE)(E. - E*)dE* E dE (B13) 

The change in the population of the j n state is derived in exactly the same way with 
the exception that equation (Bl) now represents the number of electrons leaving, and 
equation (B5) represents the number entering the bound state. Hence, 


Hence, equation (BIO) becomes 


© ir =”^(»3 
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n e( n 3 S - n j) < va I (E j ~ E *» 


(B14) 


The energy transferred to and from the electron gas due to ionization-recombination 
results in an expression similar to equation (B13). The unit mechanism implicit in the 
derivation of (B13) is denoted as 6(0f e /St)j^ where the 6 means that none of the indi- 
cated integrations in equations (Bl) and (B5) are performed; that is, the averaging is not 
yet carried out. Multiplying 6(9f e /3t)j by E and averaging gives the energy loss in 
ionization; and multiplying ^(Bfg/St)^ by E and again averaging gives the energy gained 
by recombination. Hence, 



“ e X ("> ? " "Xl mv2w i (E i - E 



(B15) 


where 


/ \ ^ j 2 

- mv 2 vu T (E. - E*)\ = 8irmf— — \ f exp(-AE)a T (E. - E*)E 2 dE dE* (B16) 

,2 1 3 / \27rm/ J 1 - 1 


Inelastic-Superelastic Collisions 

It is assumed first that Ej > E^ (or k > j). The change in population of the j state 
due to inelastic and superelastic collisions proceeds as follows: The number of particles 
that leave the j th state for the k state is equal to the inelastic collision rate (since 

E j E k^’ that is ’ 




exp(-AE)E dE 


(B16) 


where the explicit form of f g was used, and the differentials were transformed as in 
equation (B2). The inverse to this process is the superelastic collision rate, equal to 


SL, ■ -a/"- exp( -* E ‘ )E * 


(B17) 


where 
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(B18) 


E* = E - (E. - E k ) 

According to equation (1902) of reference 16, the relation between cr ex and 
is given by 


ClJ k E * CT dex = i ECT ex 

Hence, summing over all k > j gives the partial result 


(B19) 



= 87rmn 
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e \2irm./ 


3/2 YV 


z* 

k>i 


co. 


“k 


exp[(E. - E k )] - n If a ex (j - k)E exp(-AE) dE 

-HEj-Ek 


(B20) 

Observing the convention that E* is the energy of the free electron before super- 
elastic collision (or after inelastic collision) and expressing the rates in terms of the ex- 
citation cross section shows easily that for E k > E^ (or k < j) the following rate equation 
holds: 



n j 


exp[A(E k - E j)l / a ex (k - j)E exp(-AE) dE 


(B21) 


The following relations are noted: 


<v 


>0 dex> ' n 


dex dT e * 8 ™ 


X \ 3 ^ 2 /'“e- XE a dex E* dE* (B22) 

Jo 


27rm 


Using equation (B19) and the energy relation equation (B18) converts equation (B22) 


into 
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)] f* exp(-AE)a ex E dE 
E f E k 

= _j exp[A(E. - E k )]<va ex > (B23) 

k 

Retention of the asterisk in equation (B23) is spurious since both v and v* are dummy 
variables that lose their distinction after averaging. 

Adding equations (B20) and (B21) then gives the result quoted in equation (11): 


/ . \3/2 co. 


[”k - 3 s n i eXp[>(E k ' E i»} <v<, ex< k - i» 

k k<j 


+ ne X/£ ^ exptME) ’ Ek>i ’ "I 

k>i 


^< v % x (j-*k)> (B24) 


The energy exchange due to inelastic-superelastic collisions is simply computed by 
considering that for every electron excited from the state j to the state k a free elec- 
tron originally having energy E leaves the phase cell characterized by E. Thus 



= n. 



co. 


CO, 


n k exp[A(E. - E k )] 


- d H 


mv 2v a e x(j - k) 


(B25) 


j k>j 
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